A  NUMERICAL  METHOD  OF  1000-MB  PROGNOSIS 
THROUGH  PREDICTION  OF  THICKNESS  FIELD 


HO-CKIH  CHIN 

AND 
HA  SONG  PAK 

N 

us   r  UATE  SCHOOL 

UA 


A  NUMERICAL  METHOD  OF  1000-MB  PROGNOSIS 

THROUGH 
PREDICTION  OF  THICKNESS  FIELD 

Ho-Chih  Ch-  r 
Chuang-Tao  Kuan 
Ha  Song  Pak 


A  NUMERICAL  METHOD  OF  1000-MB  PROGNOSIS 

THROUGH 
PREDICTION  OF  THICKNESS  FIELD 


by 


Ho-Chih  Chin       Commander 

Republic  of  China  Navy 

Chuang-Tao  Kuan     Lieutenant  Commander 

Republic  of  China  Navy 

Ha  Song  Pak         Lieutenant  (Junior  Graded 

Republic  of  Korea  Navy 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 
METEOROLOGY 


United  States  Naval  Postgraduate  School 
Monterey,  California 

19  6  1 


'IOOL 

A  NUMERICAL  METHOD  OF  1000-MB  PROGNOSIS 

THROUGH 
PREDICTION  OF  THICKNESS  FIELD 

by 

Ho-Chih  Chin 

Chuang-Tao  Kuan 

Ha  Song  Pak 

This  work  is  accepted  as  fulfilling 

the  thesis  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 

IN 

METEOROLOGY 

from  the 

United  States  Naval  Postgraduate  School 


ACKNOWLEDGEMENTS 

The  writers  are  indebted  to  Professor  George  J. 
Haltiner  of  the  Department  of  Meteorology  and  Oceanography 
at  United  States  Naval  Postgraduate  School  for  suggesting 
the  topic,  for  outlining  the  way  in  which  the  study  might 
be  pursued,  and  for  the  assistance  he  rendered  them. 

They  are  also  indebted  to  Lieutenant  Commander  Edward 
M.  Carlstead  of  the  Fleet  Numerical  Weather  Facility, 
United  States  Navy,  for  the  assistance  he  gave  them  during 
the  course  of  investigation. 


11 


ABSTRACT 

A  two-level  numerical  prediction  model  incorporating 
terrain  and  non-adiabatic  warming  effects  is  tested.   An 
attempt  is  made  to  determine  the  empirical  coef f icient (s) 
for  each  of  three  terms  composing  a  prognostic  equation 
designed  to  predict  the  1000   to  500-mb  thickness  through 
the  use  of  a  digital  computer. 

The  model  is  applied  to  three  cases  in  April  for  a 
large  part  of  the  Northern  Hemisphere.   The  predicted 
positions  of  most  of  the  pressure  systems  were  reasonably- 
accurate;  however  the  model  tends  to  over-develop  pressure 
centers.   The  effects  of  terrain  and  non-adiabatic  warming 
were  in  accord  with  theoretical  expectations;  however,  the 
excessive  development  dominated  the  numerical  verification 
so  that  no  significant  improvement  resulted  from  the  in- 
clusion of  these  refinements. 
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1.   INTRODUCTION 

The  availability  of  high-speed  large-capacity  elec- 
tronic digital  computers  has  long  since  made  numerical 
weather  prediction  an  operational  reality.   Notable  success 
has  been  achieved  at  the  500-mb  level.   However,  surface 
prognoses  obtained  by  numerical  integration  techniques  on 
the  average  have  not  equalled  the  accuracy  obtained  by 
subjective  methods.   Even  with  some  relatively  complex 
models,  the  degree  of  success  has  not  corresponded  to  the 
degree  of  sophistication  of  the  models.   To  improve  the 
quality  of  surface  prognosis,  two  alternatives  seem  to  be 
feasible.   One  approach  is  to  devise  even  more  elaborate 
mathematical  models  in  order  to  take  into  account  as  many 
pertinent  weather  parameters  as  possible.   Perhaps  new 
parameters  which  are  not  just  a  recombination  of  the 
familiar  physical  variables — pressure,  temperature  and 
humidity,  will  have  to  be  introduced  and  measured  to  a  high 
degree  of  accuracy  on  a  synoptic  time  and  space  scale. 
Another  alternative  is  to  try  quasi-empirical  prediction 
models  with  a  reasonable  dynamical  basis.   In  this  paper 
a  simple  quasi-empirical  1000-mb  prognostic  model  is 
derived  and  testing  results  are  presented. 

Since  500-mb  prognoses  from  the  barotropic  model  or 
its  variations  have  been  quite  successful,  a  next  logical 
step  would  be  to  devise  a  method  to  forecast  the  thickness 
field;  and  then  through  the  predicted  thickness  arrive  at 
the  prognosis  of  some  other  level.   In  this  experiment,  a 
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technique  for  predicting  the  thickness  field  between 
1000  mb  and  500  mb  is  combined  with  results  of  the  baro- 
tropic  model  currently  in  use  by  the  Fleet  Numerical 
Weather  Facility.   The  predicted  500-mb  stream  field  at 
every  hour  is  applied  in  turn  to  generate  a  thickness 
between  levels  through  a  thickness-tendency  equation. 
After  sufficient  iterations,  a  24-hour  1000-mb  prognos  .0 
contour  map  is  obtained. 

From  the  hydrostatic  viewpoint,  pressure  is  simply  the 
weight  of  the  air  column  above  a  particular  level,  and 
local  changes  in  pressure  may  be  thought  of  as  the  inte- 
grated density  change  in  the  air  column.   The  density  change 
can  be  described  in  terms  of  the  temperature  field  which, 
in  turn,  can  be  expressed  as  a  thickness  change.   The  main 
feature  of  this  model,  the  thickness-tendency  equation, 
takes  into  account  the  temperature  field,  horizontal 
motions,  and  the  vertical  motions  induced  by  terrain  and 
surface  stress. 

Actually  there  is  no  inherent  restriction  on  the  choice 
of  the  500-mb  prediction  equation.   The  thickness-tendency 
equation  could  be  combined  with  any  500-mb  model,  though 
perhaps  requiring  some  minor  modifications  in  the  pro- 
gramming.  Obviously,  if  another  500-mb  model  is  employed 
instead  of  the  particular  barotropic  model  used  here,  the 
results  of  testing  may  be  quite  different,  and  the  empir- 
ical coefficients  may  not  have  the  same  values  as  obtained 
in  this  experiment. 
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2.   Prognostic  Model 

The  prognostic  model  is  composed  of  two  prediction 
equations.   The  first  is  the  500-mb  stream-barotropic  vor- 
ticity  equation,  the  other  is  a  thickness-tendency  equation, 

With  the  twisting  term  and  vertical  advection  of  vor- 
ticity  omitted,  the  vorticity  equation  in  (x,y,p,t)  coor- 
dinates becomes 

j£-  +v.v(M)=  -CM)w ,  (1) 

Due  to  the  fact  that  at  intermediate  levels  in  the  tropo- 
sphere between  500  and  600  mb,  the  horizontal  velocity 
divergence  is  normally  relatively  small  and  frequently  may 
be  neglected,  the  vorticity  equation  may  be  written  in  the 
form 


§f-  +V-7(f  +  /J=     0 


at 

Mere  geostrophic  approximation  of  (2)  with 
yields  the  barotropic  quasi-geostrophic  model 


(2) 


(3 


where  J  is  the  Jacobian  operator.  Equation  (3)  has  had 
considerable  success  in  predicting  the  height  field  for 
500  mb;  but  due  to  the  velocity  divergence  inherent  in  the 
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geostrophic  wind,  an  error  commonly  referred  to  as  "spuri- 
ous anticyclogenesis"  is  introduced.   To  remove  this  source 
of  error,  a  non-divergent  wind  is  used  as  follows: 

With  these  substitutions,  equation  (2)  becomes: 

Vz(f£)+70,  v*f+f)  =  o       .  (4) 

A  characteristic  of  the  barotropic  model  as  repre- 
sented by  (4)  is  that  very  long  waves,  say  number  1,  2,  and 
3  with  respect  to  a  latitude  circle,  tend  to  retrogress  at 
very  high  speeds,  a  phenomenon  not  observed  in  nature.   To 
control  this  error,  an  additional  term  has  been  introduced. 
The  result  is  the  Helmholtz  stream-barotropic  equation 

where  Y\       is  the  absolute  vorticity  and  JUL     is  a  constant. 

The  first  law  of  thermodynamics  and  the  Poisson's 
equation  can  be  combined  in  the  form 


Q-  ^irf It  +  vve +«. Jf )  _ 


(6) 


Also  from  Poisson's  equation  and  the  equation  of  state,  it 
follows  that 


T^Sy/cT .    ..*/<>  (v) 


0   =    dR   "p       loco 


which  may  be  used  in  (6)  to  give 


(8) 


Finally  replacing  (j^      in  the  first  two  terms  on  the  left 
by  means  of  the  hydrostatic  equation,  and  on  the  right 
through  the  equation  of  state  yields 


dt  la?/ 


-  V  V  ( 


Jf)-^ 


_  an 


(9) 


Here  fy  =™-r —  -. —  --— -  is  the  stability  factor  and  is 

w        8    trp  dp 

considered  to  be  constant  in  space  and  time. 

The  "vertical  velocity"  may  be  considered  to  consist 
of  three  contributions,  the  large-scale  vertical  motion, 
the  terrain-induced  vertical  motion  and  frictionally-in- 
duced  vertical  motion.   Here  these  contributions  will  be 
represented  by  assuming   6J   to  be  of  the  form 


to  =  "« 


^(■^fl-^^T-^VHt+wj 


(io; 


where 


1= 


rt^L   $"    ?^?^r0 


1  +  J^L  kt^rq 


y 


The  first  term  on  the  right  side  of  (10)  is  the  para- 
bolic distribution  frequently  assumed  to  be  typical  of 
large-scale  pressure  systems  {1]  . 

The  second  term  of  U)  represents  the  terrain-induced 
vertical  motion.  Here  the  assumption  is  made  that  the  sur- 
face vertical  velocity  w„  is  proportional  to  the  horizontal 
wind  and  the  gradient  of  terrain  height  at  the  surface; 
i.e.  ,  Wq  =  k^T  •  y  H..  Since  £J0  =  -g  0  wQ ,  it  follows  that 
(A-  -g  P  k  y  „^H  H..  The  terrain-generated  vertical  velo- 
city is  further  assumed  to  decrease  linearly  with  pressure; 
thus  we  arrive  at  the  form  in  the  second  term. 

The  third  term  represents  vertical  motion  induced  by 
surface  friction.   Here  CJ*  is  the  "vertical  velocity"  at 
top  of  the  friction  layer  due  to  surface  stress.   It  is 
assumed  to  have  a  simple  distribution  with  respect  to 
height  which  can  be  approximated  linearly  as  shown  in  fig. 1. 


C0=    O     o±   IP 


Cd^Qlat  f£ 


cu«  c    <&  7j 


Figure  1.   Profile  of  vertical  velocity  due  to 
surface  stress. 


Cressman  [2]  deduced  the  following  form  for  CJ^  in  terms 
of  geostrophic  surface  wind  and  a  drag  coefficient: 


to, 


i 


where 


F  = 


(12) 


Here  u  and  v  are  the  x-  and  y-components  of  the  geostrophic 
wind  at  1000  mb,  f  is  the  coriolis  parameter  and  C.  is  the 
drag  coefficient.   The  latter  is  mainly  a  function  of  the 
roughness  of  the  underlying  surface  and  stability  in  the 
micrometeorological  layer,  and  has  been  determined  for 
average  conditions  by  Cressman  for  large-scale  forecasting 
over  the  Northern  Hemisphere. 

Upon  integration  between  ~fr   and   P_   equation  (9) 


Vs>Vh  -f  it T?H  -  T%<r$(M- VHt 


becomes 

ah  _ 
at 

-Vs>vl 

— 

<r9F  ■' 
4J    i 

e*+(T,-Ti)jL. 


+  QCT1 


R 


frQ 


(13) 


where  Q  represents  the  amount  of  heat  added  to  the  layer 
per  unit  time. 

If  in  equation  (13)  we 

1.   assume  (T^£jj== 'xV'' \7  H    and  combine  it:  ™ith 
the  first  term  and  affix  a  coefficient  K-,  ; 


2.  combine  coefficients  in  the  second  term  and 
designate  the  combination  as  K- ; 

3.  assume  Q  is  proportional  to  V o0 VT  •    or 

Q  =  K-jV  «V7T  ,  where  T   is  the  sea-surface 
temperature;  and 

4.  assume  a  reasonable  pressure  level  as  top 
of  friction  layer  in  the  frictional  term 
and  designate  the  combined  coefficient  as 


V 


we  arrive  at  : 


Ik. 
at 


■=-KV5Mi-)(iM'VHt+WVVK*F    <14 


This  is  the  thickness-tendency  equation,  with  the  empirical 
coefficients  to  be  determined  through  testing.   The  thick- 
ness equation  together  with  the  500-mb  stream-barotropic 
equation  constitute  the  prognostic  model. 

The  local  rate  of  thickness  change  is  thus  seen  to  be 
due  to  the  combined  effect  of  thickness  advection  between 
layers,  terrain  effects,  non-adiabatic  heating  and  a 
frictional  effect.   Since  cold  advection  is  usually 
associated  with  adiabatic  warming  due  to  subsidence  aloft 
and  by  heating  from  below  in  the  surface  layers,  while 
warm  advection  is  usually  counteracted  by  cooling  due  to 
large-scale  upward  motion  and  in  the  surface  layers  by 

cooling  from  below,  the  coefficient  K.^  might  be  expected 
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to  have  different  values  for  cold  and  warm  advection.   This 
was  borne  out  by  the  results  of  this  numerical  experiment. 

The  terrain  term  represents  vertical  motion  induced 
by  uneven  terrain,  and  obviously  it  vanishes  over  oceans 
and  level  land. 

The  non-adiabatic  term  in  the  form  K^V  ,VT   represents 
vertical  diffusion  of  heat  between  air  and  surface  water. 
For  an  air  mass  moving  over  warmer  water,  the  surface 
temperature  quickly  adjusts  to  that  of  the  surface  water 
and  rapid  vertical  diffusion  of  heat  takes  place.   The  rate 
of  diffusion  will  normally  increase  with  increasing  hori- 
zontal wind  which  enhances  the  turbulence,  and  upward  dif- 
fusion of  heat  will  certainly  continue  as  long  as  the  water 
temperature  increases  downwind.   For  warm  air  moving  over 
colder  water  the  stability  created  by  surface  cooling  in 
the  lowest  layer  will  tend  to  inhibit  diffusive  heat  ex- 
change.  This  suggests  the  empirical  coefficient  K-.  may 
also  be  assigned  two  values  depending  on  the  sign  of 

The  effects  of  surface  stress  about  closed  pressure 
systems  may  be  summarized  as  follows: 

Low:   upward  motion  at  top  of  friction  layer,  £0  § 

negative;  contribution  to  thickness  tendency, 

negative. 
High:  downward  motion  at  top  of  friction  layer  U) fi 

positive;  contribution  to  thickness  tendency, 

positive. 


Hesselberg  and  Friedmann  found  that  the  friction  effect 
is  at  least  two  orders  of  magnitude  smaller  than  some  of  the 
other  terms  in  the  vorticity  equation.   The  assumption  is 
made  that  it  will  also  be  quite  small  as  compared  to  other 
terms  in  the  thickness-tendency  equation.   Hence  in  the 
present  experiment  it  was  neglected.   Thus  equation  (14) 
becomes 


dh 


ff  =  -KM-VK-K.Vt -VHt  +  K3V'VT,      ,», 


where  V   represents  the  wind  at  the  terrain  height.   Next 
the  winds  V,-  and  V   in  equation  (15)  are  approximated  by 
geostrophic  values  resulting  in 


&=  -  jr  [KiT(*,.«-KlT(*.;i)]-  KJfrVHt™ 

The  wind  at  terrain  height  is  computed  by  linear  inter- 
polation.  On  the  assumption  that  wind  hodograph  is  approx- 
imately a  straight  line  between  1000  and  500  mb ,  V.  may  be 
expressed  in  the  form 


vt=M-(^ 


Ik 


T     .  (17) 
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Thus  the  wind  at  terrain  height  is  equal  to  the  wind  at 
500-mb  level  minus  some  fraction  of  the  thermal  wind,  and 
the  fraction  is  the  ratio  of  the  surface/500-mb  thickness 
to  the  500/1000-mb  thickness.   With  the  geostrophic  approx- 
imation for  the  wind  and  the  letter  A  to  designate  ( — ^r~       ) 
it  follows  that 


M- 


!Y    ^+#jv 


5  \    =»1 


jx 


-Af  (-J^  +  * 


(18) 


Hence  the  terrain  term  in  equation  (16)  may  be  computed 
using  two  Jacobians.   However,  because  of  limited  computer 
storage  capacity, instead  of  the  two  Jacobian  operations,  the 
following  transformation  is  made 


V^H-f-fef  ^  lfi)-A(-|i  f  M-VHt 


(19) 


In  finite  differance  form,  the  last  expression  becomes 


Vt;VHt=%^a-(^-A^)+(2«-Ah^HtE-HtN) 


+ 


'( ^E-A  M  "  ( %s  w~  A  h  Jl  (Ht  N-Ht5)l 


)   (20) 
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where  the  subscripts  N,  S,  E  and  W  are  notations  for  the 
location  of  grid  points.   The  quantity  A  has  three  var- 
iables, but  may  be  simplified  by  taking  standard  atmos- 
pheric values  for  Z^   and  h  without  introducing  too  large 
an  error.   With  the  letters  ZfcN,  Z.-,  Z    and  Z,w  to  denote 
the  quantities  within  the  brackets  in  the  final  expression 
for  Vt*VH4.  above,  equation  (15)  has  the  form 

jf  =  ~  ^le  I K  y  fa,  h.) "  K  J  (2,  Js) 


+K  Jtetj-^XHu-  HtwHC^u  -2  J(HtfrHt 
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3.   Procedure 

Three  days'  data  in  April  1955  were  tested.   Progno- 
ses were  made  for  April  2nd,  3rd  and  4th  and  were  labeled 
as  Cases  1,  2  and  3,  respectively   A  CDC-1604  digital 
computer  was  employed  for  the  computations  over  a  197  7 
point  octagonal  grid  for  the  Northern  Hemisphere  poleward 
of  ION.   At  the  end  of  each  iteration,  a  predicted  thick- 
ness is  obtained  at  each  grid  point.   This  is  subtracted 
from  the  prognostic  500-mb  stream  function  to  find  the 
1000-mb  height.   The  final  output  is  a  predicted  1000-mb 
height  field  in  a  contour  map  form.   The  forward-difference 
method  was  used  for  calculation  of  ~     for  the  first  hour 
and  a  centered-dif ference  method  was  used  for  succeeding 
hours.   The  forward  and  centered  difference  methods  may  be 
expressed  respectively  as 


i*+.t-  ''■'Hiil** 


and 


hf-hAt—      ^t-Af  +   (J£Tjt^   A   *" 


A  schematic  block  diagram  of  the  thickness  prediction 
program  is  given  in  fig.  2. 

Using  the  computer,  a  "pillow"  and  a  RMSE  for  each 
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Figure  2.   Schematic  block  diagram  showing  procedure, 
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prognostic  map  were  computed  by  comparison  with  the  actual 
map  for  verification  purposes.   The  formulae  used  in  com- 
puting the  pillow  and  RMSE  are 


pillow   = 


Yl-  I 


(A-B) 


x 


RMSE     =    -Ja^1- 


M-B)-^lUui]^ 


X 


where  A  and  B  stand  for  predicted  and  verifying  values, 
respectively,  at  the  same  grid  point;  and  x  stands  for  the 
total  number  of  points  differenced.   Qualitatively,  a  pillow 
is  the  algebraic  mean  difference  over  the  whole  grid  between 
the  two  fields  A  and  B?  while  the  RMSE  is  a  numerical  meas- 
ure of  the  prognosis,  and  it  reflects  errors  in  position 
and  in  intensity  of  pressure  systems. 

Due  to  fictitious  reflections  on  the  boundary  points 
and  inaccuracies  in  the  geostrophic  approximation  at  low 
latitudes,  it  was  decided  to  limit  the  verification  to  the 
areas  north  of  the  20th  latitude  circle. 

Since  the  major  contribution  to  thickness  change  is 
due  to  the  advection  term,  it  was  investigated  first. 
Values  of  K,  determined  by  previous  investigators  [3]  were 
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used  initially,  then  the  value  was  varied  to  arrive  at  the 
best  verification.   Furthermore,  since  warm  and  cold  ad- 
vection  did  not  necessarily  give  the  same  contribution 
to  the  predicted  change  in  thickness,  combinations  of  two 
different  values  were  tried  during  the  investigation  of 
this  term. 

With  the  coefficients  for  the  advection  term  tenta- 
tively determined,  the  investigation  was  carried  on  to 
find  suitable  coefficients  for  the  terrain  term.   Only  one 
coefficient  of  best  fit  was  sought  for  the  terrain  term. 
During  the  investigation  of  this  term,  a  smoothed  terrain- 
height  field  is  required  in  order  to  obtain  the  signifi- 
cant scale  of  vertical  motion.   Since  the  grid  distance  is 
about  380  km,  any  smaller-scale  variation  in  the  terrain 
field  would  not  be  included.   For  this  reason,  the  smoothed 
terrain-height  field  data  of  J.  Smagorinsky  was  used. 

Finally,  an  investigation  was  conducted  to  find  suit- 
able coefficients  for  the  non-adiabatic  term.   Two  coef- 
ficients were  sought  for  this  term,  one  for  heating  and 
one  for  cooling.   Values  of  K  determined  by  previous  inves- 
tigators [4]  were  also  used  initially  in  this  case.   The 
monthly-mean  sea-temperature  field  of  April  was  used  during 
the  investigation  of  the  non-adiabatic  term.   This  field 
was  obtained  by  graphical  addition  of  the  monthly-mean 
air-temperature  charts  of  April  and  the  3-month-mean  air- 
sea  temperature  difference  charts  [5] .   The  resulting 
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temperature  field  is  shown  in  fig.  10.   The  temperature 
field  on  the  land  area  was  included  to  facilitate  the 
computations;  however,  these  values  did  not  enter  into  the 
final  results. 
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4.   Results  and  Discussion: 

Figures  3  to  8  indicate  the  initial  maps  and  the 
prognostic  maps  using  different  combinations  of  the  terms 
forming  the  thickness-tendency  equation.   The  isolines 
are  drawn  at  intervals  of  200  feet,  and  labeled  in  hun- 
dreds of  feet.   Case  3  has  been  chosen  for  illustration 
purposes  because  it  gave  the  best  results  among  the  three 
cases  tested,  both  in  the  prediction  of  the  positions  of 
the  pressure  centers  and  in  numerical  verifications. 

By  comparing  figs.  4  and  5,  it  can  be  seen  that  the 
positions  of  the  pressure  systems  are  in  close  agreement 
with  four  exceptions,  namely:   the  high  over  the  east  coast 
of  Asia;  the  high  off  the  west  coast  of  the  United  States; 
the  high  over  the  Caspian  Sea;  and  the  high  lying  to  the 
west  of  Gibraltar.   Possible  reasons  for  errors  in  their 
locations  will  be  discussed  at  the  end  of  this  section. 
It  was  also  of  interest  to  notice  that  the  positions  of 
systems  in  the  Atlantic  Ocean  were  predicted  more  accurately 
than  those  in  the  Pacific  Ocean. 

Before  discussing  the  coefficients  found  thus  far  in 
this  investigation,  it  is  necessary  to  mention  two  phenom- 
ena related  to  development.   These  phenomena  are  fictitious 
anticyclogenesis  and  cyclogenesis  which  predominated  in 
all  three  cases  tested.   While  no  dynamical  explanation 
could  be  determined  from  this  investigation,  it  was  apparent 
that  the  effect  was  due  mainly  to  the  advection  term. 
Comparison  of  figs.  4  and  6  shows  that  most  of  the  highs 
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were  too  high  by  an  average  amount  of  200  feet,  while  most 
of  the  lows  were  too  low  by  an  average  amount  of  100  feet 
when  the  advection  term  was  used  alone  to  make  the  forecast. 
Since  the  major  contribution  to  thickness  change  is  due  to 
the  advection  term,  these  phenomena  invariably  concealed 
the  effects  of  a  refinement  like  the  non-adiabatic  heat 
exchange  in  the  numerical  verification. 

Aside  from  the  effect  of  over-prediction  by  the  advec- 
tion term,  it  was  found  that  the  combination  of  coeffi- 
cients 0.4  for  warm  advection  and  0.5  for  cold  advection 
gave  the  best  verification  when  this  term  is  used  alone  to 
make  the  forecast.   Other  sets  of  coefficients,  such  as 
0.8/0.9,  0.7/0.8,  0.5/0.5,  0.3/0.4,  to  mention  only  a  few, 
gave  larger  RMSE  values.   Therefore,  it  was  decided  tenta- 
tively to  use  0.4/0.5  as  the  coefficients  of  best  fit  for 
the  advection  term.   The  results  of  numerical  verification 
for  prognoses  using  the  advection  term  only  with  coeffi- 
cients of  0.4  for  warm  advection  and  0.5  for  cold  advection 
are  shown  in  Table  1. 
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Case 

Pillow 

RMSE 

Number 

(ft) 

(ft) 

1 

(OOOOZ  April  2, 

0 

-f  189 

1955       ) 

2 

(OOOOZ  April  3, 

+  14 

4-  171 

1955        ) 

nm^r^T,- 

3 

(OOOOZ  April  4, 

+  7 

-|~168 

1955        ) 

Table  1.   Results  of  numerical  verification  for  prognoses 
using  the  thickness  advection  term  only  with 
coefficients  of  0.4  for  warm  advection  and  0,5 
for  cold  advection. 
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For  the  investigation  of  the  terrain  term,  verifica- 
tion by  pillow  and  RMSE  in  Table  2  showed  no  improvement 
in  RMSE  when  the  terrain  term  was  included  in  this  model. 
In  this  connection,  it  should  be  noted  that  the  empirically- 
determined  best-fit  coefficient  of  the  terrain  term  is  0.01, 
while  the  maximum  terrain-height  gradient  is  about  12  times 
larger  than  the  thickness  gradient.   Hence  the  terrain  term 
contribution  can  be  at  most  about  25%  of  the  advection  term. 
Also,  less  than  half  of  the  total  area  covered  by  the  grid 
point  is  affected  by  the  terrain  term,  while  the  pillow 
and  RMSE  were  computed  over  the  entire  grid.   Therefore, 
the  actual  contribution  of  the  terrain  term  should  be 
larger  than  that  shown  by  RMSE. 

In  order  to  investigate  the  effects  of  the  terrain- 
induced  vertical  motion  in  more  detail   a  differential 
analysis  was  made  and  presented  in  fig.  9.   Here  the  dif- 
ference between  the  forecast  with  thickness  advection  only 
(fig.  6)  and  that  with  both  the  terrain  term  and  thickness 
advection  (fig.  7)  was  computed.   The  isolines  are  drawn 
at  intervals  of  50  feet,  and  the  figures  labeled  in 
tens  of  feet.   Of  the  20  closed  isolines  observed  in  fig.  9, 
14  indicate  a  decreased  error;  5  of  them  an  increased  er- 
ror; and  one,  neutral.   Thus,  there  are  some  areas  where 
the  terrain  term  does  not  give  any  improvement  and,  in 
fact,  made  the  results  worse.   For  example,  over  Greenland, 
the  forecast  made  with  the  thickness  advection  alone  over- 
intensifies  the  low,  and  with  the  inclusion  of  the  terrain 
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Case 

Number 

Pillow 
(ft) 

RMSE 
(ft) 

1 

(OOOOZ  April  2, 
1955       ) 

-  15 

+  191 

■ 

2 

(OOOOZ  April  3, 
1955       ) 

+  1 

t  172 

3 

(OOOOZ  April  4, 
1955       ) 

-  7 

+  168 

Table  2.   Results  of  numerical  verification  for  prog- 
noses using  advection  and  terrain  terms  with 
coefficients  of  0.4  for  warm  advection,  0.5 
for  cold  advection  and  0.01  for  the  terrain 
term. 
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term,  further  deepening  is  observed  on  the  leeward  side  of 
the  mountain,  giving  an  even  poorer  result. 

Fourteen  closed  isolines  of  the  20  on  the  differential 
analysis  showed  a  better  forecast,  with  a  two- to  eight-mb 
improvement  in  the  surface  pressure  field,  when  the 
terrain-induced  vertical  motion  over  land  areas  was  in- 
cluded.  These  results  indicate  some  justification  for  the 
inclusion  of  the  terrain  term  in  the  model. 

Since  the  terrain-induced  vertical  motion  is  a  func- 
tion of  the  wind  component  parallel  to  the  terrain-height 
gradient  as  well  as  the  terrain-height  gradient  itself, 
large  values  of  terrain  contribution  may  be  expected  in 
areas  where  above-mentioned  values  are  large.   In  this 
respect,  the  test  data  available  for  this  investigation 
were  not  favorable  for  a  good  evaluation  of  the  terrain 
term.   The  synoptic  maps  show  that  the  areas  where  the 
terrain  gradient  is  large  are  mainly  associated  with 
small  wind  components  parallel  to  the  terrain-height 
gradient. 

With  only  three  days'  data  tested,  no  final  conclu- 
sion can  be  reached  regarding  the  terrain-induced  verti- 
cal motion  in  this  model  at  this  time.   Effects  of  the  ter- 
rain-induced vertical  motion  were  produced  in  areas  where 
they  were  expected  from  dynamic  considerations;  however,  the 
values  are  too  small  in  general.   This  seems  to  suggest 
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that  further  experiments  may  indicate  a  larger  coeffi- 
cient for  the  terrain  term  with  perhaps  some  other  addi- 
tive empirical  factors. 

The  results  of  numerical  verification  for  prognoses 
using  advection  and  non-adiabatic  terms  are  shown  in  Table 
3.  The  eight  blanks  in  Table  3  arose  in  cases  where  the 
forecast  map  was  very  much  distorted,  in  contrast  to  the 
regular  appearance  of  the  remaining  forecast  maps,  This 
erratic  behavior  was  apparently  connected  with  a  malfunc- 
tion of  the  magnetic  tape  units. 

The  coefficients  of  best  fit  for  the  non-adiabatic 
term,  as  shown  in  Table  4,  were  10  for  heating  and  6  for 
cooling.   They  were  chosen  only  because  they  provided  the 
smallest  pillow.   The  RMSE  obtained  by  using  this  set  of 
coefficients  was  -j-  168  feet  for  case  3,  which  equaled 
exactly  that  obtained  by  using  advection  term  alone  for 
prognosis.   As  a  matter  of  fact,  Table  3  shows  that  the 
RMSE  values  increased  when  the  coefficients  were  increased 
algebraically.   Since  it  was  mentioned  earlier  that  the 
RMSE  is  a  numerical  measure  of  the  accuracy  of  the  forecast, 
it  appears  that  the  addition  of  the  non-adiabatic  term  at 
best  made  no  improvement  in  the  forecast.   Two  arguments 
need  to  be  made  clear  at  this  point.   Firstly,  fictitious 
anticyclogenesis  and  cyclogenesis  by  the  advection  term 
alone  already  have  over-intensified  most  of  the  highs  and 
the  lows.   Since  the  effect  of  the  non-adiabatic  term  on 
the  average  is  one  order  of  magnitude  less  than  that  of 
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Coeff . 

Case 

1 

Case 

2            Case  3 

heat- /cool- 

(0000Z Apr 

2.   '55) 

(0000Z  Apr 

3, '55) 

(0000Z  Api 

c  4,  '  55) 

ing  /  ing 

Pillow 

RMSE 

Pillow 

RMSE 

Pillow 

RMSE 

/ 

(ft) 

- — _., ,  . 

(ft) 

(ft) 

(ft) 

(ft) 

(ft) 

5/5 

-15 

+  192 

+  1 

-t  175 

-  4 

+  167 

5/10 

-  6 

+  200 

+  9 

•\   178 

+  3 

+  170 

5/15 

-  4 

4  205 

+  17 

+  182 

+  12 

+  178 

5/20 

■f  9 

+  209 

+■  24 

+  188 

+  20 

+  187 

10/5 

-  4 

+  199 

-  4 

+  177 

-11 

1 
t  168 

10/6 

-  2 

-f-200 

-  3 

+  177 

-  1 

+  168 

10/10 
10/15 
10/20 

-  6 

+  215 

f  10 
+■  10 
+  18 

+  183 
+  184 
4  189 

+  4 
+  5 

+  13 

1 

+  177 
+  178 
4  187 

0 

\-  219 

15/5 
15/10 

-11 
-  3 

i 

+■  182 
+  184 

-18 
-  9 

+  171   | 
+  175 

-23 

+  218 

15/15 

-13 

-1-  232 

+  4 

+  188 

o 

1 

+  182 

15/20 

-  8 

+  242 

+  11 

+  193 

+  7 

t  192  1 
] 

20/5 

-43 

f  227 

-18 

+  190 



20/10 

-32 

+  235 

-10 

+  191 

-15 

+  181 

20/15 

-22 

+  248 

"  2 

+  196 

-  7 

*■  190  ! 

20/20 

+  5 

+  200 

+  1 

+  202 

25/5 

-42 

i 

f  230 

• 

-20 

+  191 

> 

|   -24 

+  179 

25/10 

-33 

t  238 

1 

-10 

+  193 

-16 

+  182 

25/15 

-23 

*  250 

"  3 

+  197 

i 

-  8 

i 

+  192 

25/20 

-13 

+  270 

+  4 

■w  201 

0 

+  204 

(Continued  on  following  page) 
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(Continued  from  previous  page) 


Coeff 

heat- /cool- 
ing    /     ing 


Case    1 
(0000Z   Apr    2, '55) 


Pillow 
(ft) 


RMSE 
(ft) 


Case  2 
(0000Z  Apr  3, '55) 


Pillow 
(ft) 


RMSE 
(ft) 


Case  3 
(0000Z  Apr  4.,  '  55) 


Pillow 
(ft) 


RMSE 
(ft) 


30/5 
30/10 
30/15 
30/20 


■40 
■24 
•16 
•  8 


-i-  218 
+  219 
r  221 
r  228 


43 

-t-  200 

30 

-f-  204 

20 

1-  218 

12 

t  226 

Table  3.   Results  of  numerical  verification  for  prognoses 
using  advection  and  non-adiabatic  terms  with 
coefficients  of  0.4  for  warm  advection  and  0.5 
for  cold  advection.   The  coefficients  used  for 
the  non-adiabatic  term  are  shown  in  the  table. 
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Case 
Number 


1 
(OOOOZ 
April  2,  1955) 


Coefficients 
used 


K 


(Warm       0.4 
advection) 

1(Cold        0.5 
advection) 


K. 


0.01 


K  (Heating)   10 
3 (Cooling)    6 


Pillow 
(ft) 


RMSE 
(ft) 


-  22 


+  199 


(OOOOZ 
April  3.  1955) 


K 


(Warm        0.4 
advection) 

1(Cold        0.5 
advection) 

K„  0.01 


K 


(Heating)   10 
3 (Cooling)    6 


-  10 


+  178 


3 
(OOOOZ 
April  4,  1955) 


Ki 


K, 


(Warm       0.4 
advection) 

(Cold       0.5 
advection) 


0.01 


K. 


(Heating)   10 
(Cooling)    6 


-  9 


-I-  168 


Table  4.   Results  of  numerical  verification  for  prog- 
noses using  all  three  terms  of  the  thickness- 
tendency  equation,  with  coefficients  as  shown 
in  the  table. 
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the  advection  term,  the  advantage  of  including  the  non- 
adiabatic  term  would  hardly  appear.   Secondly,  for  the  non- 
adiabatic  heat  exchange,  the  zones  of  thermal  concentration, 
where  transformation  of  polar-continental  air  to  polar- 
maritime  air  takes  place,  are  really  of  primary  concern. 
During  the  month  of  April,  these  zones  of  thermal  concentra- 
tion were  found  along  the  northeast  coasts  of  Asia  and  North 
America.   For  the  three  cases  investigated,  cyclonic  activ- 
ity prevailed  in  these  areas.   With  the  1000-mb  lows  al- 
ready over-intensified  by  the  advection  term  and  with  the 
flow  oriented  such  that  heating  took  place  in  these  areas, 
the  lows  were  intensified  further.   For  example,  the  low 
off  the  east  coast  of  Japan  in  case  3,  where  the  lowest 
observed  height  was  +■  30  feet,  the  prognosis  with  the  ad- 
vection term  alone  gave  a  lowest  height  of  -  170  feet  while 
the  non-adiabatic  term  decreased  the  lowest  value  further 
to  -  250  feet.   Thus,  in  spite  of  the  fact  that  the  behav- 
ior was  in  accord  with  theoretical  expectations,  the  RMSE 
did  not  improve  when  the  non-adiabatic  term  was  included 
in  this  case. 

In  order  to  investigate  the  effect  of  the  non-adiabatic 
term  further,  a  differential  analysis  was  made  by  sub- 
tracting a  prognostic  map  with  both  advection  and  non- 
adiabatic  terms  included  (fig.  8)  from  the  one  with 
only  the  advection  term.   This  isolated  the  effect  of  the 
non-adiabatic  term  (fig.  6)  and  it  was  found  that  by  using 
a  coefficient  of  10  for  heating,  the  largest  24-hr 
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contributions  by  this  term  in  areas  of  thermal  concentra- 
tion were  -282  feet  off  the  east  coast  of  Nova  Scotia,  at 
41N,  55W;  and  -230  feet  off  the  west  coast  of  Japan,  at  40N 
134W. 

Two  possible  sources  of  error  warrant  discussion.   In 
the  first  place,  the  1000-mb  map  was  a  relatively  flat  map, 
with  the  range  of  thickness  values  and  500-mb  height  much 
greater  than  the  range  of  the  1000-mb  height.   Therefore, 
errors  in  the  500-mb  height,  when  transmitted  downward 
through  the  thickness,  can  cause  serious  deviations  at  the 
lower  surface.   Furthermore,  it  was  obvious  that  one  could 
not  expect  the  locations  of  the  pressure  systems  on  the 
lower  surface  to  be  better  than  what  the  500-mb  prognosis 
can  provide  for  this  model.   Secondly,  the  fact  that  the 
positions  of  the  pressure  systems  were  all  better  in  the 
Atlantic  Ocean  than  those  in  the  Pacific  Ocean,  true  in  all 
three  cases,  was  most  probably  due  to  the  fact  that  weather 
reports  in  the  Atlantic  Ocean  were  more  numerous  and  accu- 
rate than  those  in  the  Pacific  Ocean.   In  other  words,  the 
initial  data  from  which  the  prognosis  was  started  could 
have  been  in  error  to  a  certain  degree  in  the  Pacific  Ocean. 
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5.   Conclusions: 

Starting  with  a  dynamical  basis  and  utilizing  the 
rapid  computing  capabilities  of  the  CDC-1604  computer, 
empirical  coefficients  were  found  for  each  of  the  three 
terms  of  the  prognostic  eguation  designed  to  predict  the 

« 

1000/500-mb  thickness.   By  subtracting  the  predicted 
thickness  from  the  500-mb  24-hr  prognosticated  stream 
function,  a  24-hr  prognosis  of  the  1000-mb  map  was  obtained, 
The  fact  that  this  model  did  provide  a  map  with  most  of  the 
pressure  systems  appearing  at  the  right  places  is  encourag- 
ing.  Due  to  the  relatively  small  sample  tested,  the  empir- 
ical coefficients  found  thus  far  were  not  necessarily  the 
best  ones;  but  they  will  nevertheless  provide  a  starting 
point  for  further  investigations  of  a  similar  nature. 

Improvements  of  the  model  might  include  the  incorpor- 
ation of  certain  empirical  rules  in  order  to  remedy  the 
phenomena  of  over-predicting  cyclogensis  and  anticyclo- 
genesis  and  the  use  of  the  weekly-mean  sea- temperature 
field  instead  of  the  monthly  mean.   Furthermore,  any  im- 
provement in  the  prediction  of  the  500-mb  maps  will  give 
increased  accuracy  at  1000  mb„   In  conclusion,  it  is  the 
opinion  of  the  present  investigators  that  this  model  can 
provide  a  first  approximation  in  the  preparation  of  the 
surface  prognostic  map. 
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Tigure    3.  The   1000-mb  map   for   0000Z  April    3,    1955 
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Figure  4.     The  observed  map  for  0000Z  April  4,  1955 
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Figure  5.     The  24-hr  prognostic  map  with  all  three 
terms  of  the  prognostic  equation  in- 
cluded. 


32 


Figure  6. 


The  24-hr  prognostic  map  with  the  thickness 
advection  term  only. 
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Figure  7.     The  24-hr  prognostic  map  with  the 

thickness  advection  and  the  terrain 
terms. 
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Figure  8. 


The   24-hr  prognostic  map  with  the  thick- 
ness advection  and  the  non-adiabatic 
terms . 
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Pigure  9.     Analysis  of  the  effects  of  the  terrair 
term. 
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Figure  10.   The  monthly-mean  sea-temperature  field 
of  April. 
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